SOBOLEV REGULARITY FOR MONGE- AMPERE TYPE EQUATIONS 



GUIDO DE PHILIPPIS AND ALESSIO FIGALLI 

Abstract. In this note we prove that, if the cost function satisfies some necessary structural 
conditions and the densities are bounded away from zero and infinity, then strictly c-convex po- 
tentials arising in optimal transportation belong to W^ +K for some k > 0. This generalizes some 
recents results [9l 1101 [22] concerning the regularity of strictly convex Alexandrov solutions of the 
Monge- Ampere equation with right hand side bounded away from zero and infinity. 



1. Introduction 

Let f2 C M. n be a bounded open set. We want to investigate the regularity of solutions to 
Monge- Ampere type equations of the form 

(1.1) det (D 2 u - A(x, Du)) = f in U, 

where / > and A(x,p) is a n x n symmetric matrix. 

This class of equations naturally arises in optimal transportation, and in reflector and refractor 
shape design problems. In these applications, the matrix A and the right and side / are given by 



A(x,Du(x)) = -D xx c(x,T u (x)), f(x) = \det(D xy c(x,T u (x)))\ 



pi(T u (x)) 



where c(x, y) represents the cost function, po and p\ are probability densities, and T u is the optimal 
transport map sending po onto p\. Under a twist assumption on the cost (see (C2) below), the 
map T u is uniquely determined through the relation 

-D x c(x,T u (x)) = Du(x). 

Moreover, when A = the above equation reduces to the classical Monge- Ampere equation. 

The regularity for the above class of equations has received a lot of attention in the last years 
[TTJ [121 El HI Q21 [ZD [231 123] - In particular, under some necessary structural conditions on A (see 
(CI) below), one can show that if / is smooth then u is smooth as well |18 [ l21 [[23ll24| . In addition, 
it is proved in [11] that solutions are locally C 1,a when / is merely bounded away from zero and 
infinity (see also (121 [T7] ) . 



Recently, the authors introduced new techniques to address the Sobolev regularity of u when 
A = 0: more precisely, under the assumption that / is bounded away from zero and infinity it 
is proved in [9] that D 2 u S LlogLi oc (f2), and with a variant of the same techniques this result 
has been improved in [10] to u £ Wi ' c +K (^) for some k > (see also [22]). Let us mention that 
these results played a crucial role in [H [2] to show the existence of distributional solutions to the 
semi-geostrophic system. 

l 
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The aim of this paper is to extend the W^' 1+K regularity to the general class of Monge- Ampere 
equations in (jl.ip . Apart from its own interest, it seems likely that this result could have applica- 
tions in the study of generalized semi-geostrophic system on Riemannian manifolds [B] , in particular 
on the sphere [20] and its perturbations EJ HU E] . 

In order to describe our result, we need to introduce some more notation and the main assump- 
tions on the cost functions. 

Let X C R n be an open set, and u : X — > R be a c-convex function, i.e., u can be written as 

(1.2) u(x) = max{— c(x, y) + Xy] 

y&Y 

for some open set Y C R n , and X y G R for all y G Y. We are going to assume that u is an 
Alexandrov solution of (jl.ip inside some open set $7 C X, i.e., 

\d c u(E)\ = [ f for all E C fl Borel, 
Je 

where 

d c u(E) := |J d c u(x), d c u(x) := {y G F : = -c(x, y) + X y }, 

xeE 

and |F| denotes the Lebesgue measure of a set F. It is well-known that, in order to prove some 
regularity results, (jl.ip needs to be coupled with some boundary conditions: for instance, when 
A = and / = 1, solutions are smooth whenever they are strictly convex, and to obtain strict 
convexity some suitable boundary conditions are needed [31 [5] . 

For the general case in (II. ip . let u be a c-convex function associated to an optimal transport 
problem, and for any y £Y define the contact set 

A y := {x G X : u(x) = —c(x, y) + X y }. 

Under some structural assumptions on the cost functions (which we shall describe below) and some 
convexity hypotheses on the supports of the source and target measure, it has been proved in 
|llj that u is an Alexandrov solution of (jl.ip inside X, and it is strictly c-convex (i.e., for any 
y G d c u(X) the contact set Ay reduces to one point) provided / is bounded away from zero and 
infinity. 

Here, since we want to investigate the interior regularity of u, instead of assuming that u comes 
from an optimal transportation problem where the supports of the source and target measure 
enjoy some global "c-convexity" property, we work assuming directly that u is a strictly c-convex 
Alexandrov solution near some point x G X, and we prove regularity of u in a neighborhood 
of x. This has the advantage of making our result more general and flexible for possible future 
applications. 

Hence, we assume that there exist (x,y) G X xY such that A^ = {x}, we consider a neighborhood 
£1 of x given by 

(1.3) n :={x e X : u(z) < -c(x, y) + X y + 5}, 

where 5 > is a small constant chosen so that Q CC X and d c u(£l) CCY (such a constant 5 exists 
because A^ := {x}). Also, we assume that u is an Alexandrov solution of 

, 14 s f det (D 2 u - A(x, Du)) = f in U, 

\ u = —c(-, y) + const on 9fl. 
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Before stating our result, let us introduce the main conditions on the 
above, and let O CC Y be a open neighborhood of d c u(Q). We define 

(1-5) |||c||| := ||c|lc3(n X 0) + \\D X xyyc\\ L oc,(Q X Q-j + ||log | det D xy c\\ 

and assume that the following hold: 

(CO) |||c||| < oo. 

(CI) For every x £ Q and p := —D x c{x,y) with y £ O, it holds 

(1.6) D PkPi Ai3(x,p)tej7] k ru>0, V£,??er\ C-t? = 0, 

where A is defined through c by Aij(x,p) := —D XiXj c(x,y), and we use the summation 
convention over repeated indices. 

Let us point out that, up to reduce the size of f2 and O (this is possible because — > {x} 
and d c u(£l) — > {y} as 5 — > 0), as a consequence of (CO) (more precisely, from the fact that 
det D xy c(x, y) ^ and by the implicit function theorem) we can assume that the following holds: 

(C2) For every (x,y) £ 0x0, the maps x £ $7 h- > —D y c(x,y) and y £ i-> —D x c(x,y) are 
diffeomorphisms on their respective ranges. 

We also notice that, because of the boundary condition u = — c(-,y) + const on 90, if / is 
bounded away from zero and infinity inside fi, then any c-convex Alexandrox solution of (II. 4h is 
strictly c-convex inside (this is an immediate consequence of [111 Remark 7.2]). Here is our result: 

Theorem 1.1. Let u : O — > M be a c-convex Alexandrov solution of (|1.4j) . Assume that c satisfies 
conditions (C0)-(C2), and that < A < / < 1/A. T/ien it G W^' c 1+ft (0) for some k>0. 

Theorem 11.11 generalizes the corresponding result for the classical Monge- Ampere equation to 
the wider class of equations considered here. With respect to the arguments in [91 [10], additional 
complications arise from the fact that, in contrast with the classical Monge- Ampere equation, in 
general (jl.ip is not affinely invariant. 

Acknowledgements: AF is partially supported by NSF Grant DMS-0969962. Both authors ac- 
knowledge the support of the ERC ADG Grant GeMeThNES. The first author thanks the hos- 
pitality of the Mathematics Department at the University of Texas at Austin, where part of this 
work has been done. 

2. Notation and preliminary results 

Through all the paper, we call universal any constant which depends only on the data, i.e., on n, 
0, O, A, and |||c|||. We use C to denote a universal constant larger than 1 whose value may change 
from line to line, and we use the notation a ~ b to indicate that the ratio a/6 is bounded from 
above and below by positive universal constants. 

An immediate consequence of the definition of c-convexity (jl.2p is that, for any xo £ X, there 
exists yo € Y such that 

u(x) > -c(x, y ) + u(x ) + c(x , y ) V x G X, 

and in this case yo £ d c u(xo). If in addition u E C 2 , then it is easily seen that Du(xq) = 
—D x c(xo,yo) and D 2 u(x$) > —D xx c(xo,yo) = A(xq, Du(xo)), where A is defined in (CI) above. 
In particular equation (|1.4p is degenerate elliptic when restricted to c-convex function. 



cost function: let il be as 



iL°°(nxe>) ' 
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It has been discovered independently in [TT] and [T7] that, because of (CI), for any xq £ f2 and 

yo G d c u(xo), through the change of variables x \-t q(x) := — D y c(x,yo) the function 

(2.1) u(q) := u(x(q)) + c(x(q), y ) - u(x ) - c(x , yo) 

has convex level sets inside f2 (here and in the sequel x(q) denotes the inverse of q(x), which is well 
defined because of (C2)). Moreover u is c-convex, where 

(2.2) c(q, y) := c(x(q), y) - c(x(q), y ), 

see [TTJ Theorem 4.3]. 

Since u solves (|1.4|) one can check by a direct computation that u solves 

(2.3) det (D 2 u - B(q, Duj) = g, 
with 

(2.4) Bij(q,p) = -D qiqj c(q,Tu(q)) and g(q) = f(x(q)) [det D xy c(x(q), T a (q))]~ 2 , 

where T s is the map uniquely identified by the relation Du(q) = —D q c(q,T a (q)). Moreover it holds 

(2.5) B lJ (-,0)=0, D p B lJ (-,0) = 0, 
so using Taylor's formula we can write 

(2.6) Bij (q, Du) = B ijM (q, Du)d k udiu 
where 

(2.7) B ijM (q, Du{q)) :=/ D PkPl Bij(q,TDu(q)) dr. 

Jo 

In addition, since condition (CI) is tensorial \21\ \19\ [T6] and |||c||| involves only mixed fourth 
derivative, it is easily seen that |||c||| ~ |||c||| and B satisfies the same assumptions as A. In particular 
(CI) and (gT]) imply that 

(2.8) B ijM CiCjVkV£>0 Ve-7? = 0. 



Given a C 1 c-convex function as above, for any xq £ fi, yo = ^(^0)1 an d h S M + , we define the 
section centered at xq of height h as 

S%(x ) := {x £ Q : u(x) < -c(x, y ) + u(x ) + c(x , y ) + h}. 

Assuming that S^(xo) C fi, through the change of variables x q(x) := —D y c(x,yo) this section 
is transformed into the convex set 

QUq Q ) := -D y c(S%(x ),y ) = {q : u{q) < h}. 

When no confusion arises, we will often abbreviate Sh(xo) and Qh(lo) for S^(xo) and Q^(qo)- 

We also recall [15] that, given an open bounded convex set Q, there exists an ellipsoid E such 
that 

(2.9) E CQ CnE, 

where the dilation is done with respect to the center of E. We refer to it as the John ellipsoid of 
Q, and we say that Q is normalized if E = B(0, 1). An immediate consequence of (|2.9|) is that 
any open bounded convex set Q admits an afline transformation L such that L(Q) is normalized. 
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Hence, given u and S^(xq) as above, we can consider u, his section Q^qo), and the normalizing 
affine transformation L. Then we define w : L(Qh) — > M as 

(2.10) w{q') := (det L) 2 / n u{q), q' := Lq. 
It is easy to check that w solves 

(2.11) det {D 2 w(q') - C(q>, Dw(q'))) = g{lT^), 
where 

C(q',Dw(q')) := (det L) 2 /™(L*) _1 £> ( L~ x q' , (det L)~ 2 / n L* Dw{q')) L' 1 



Up to an isometry, we can assume that 

n „2 



( 2 



with r\ < . . . < r n . Then L 1 = diag(ri, . . . , r n ), and 

(2.12) C ij (q\Dw(q'))=C ljM (q',Dw(q'))d k wd e w 
with 

(2.13) C ijM (q', Dw(q')) = (n . . . r n ) 2 /"^-% fc£ (<7, 

see (|2.7h . Moreover, by (12. 8p (or again because of the tensorial nature of condition (CI)) 

(2.14) C ijM ^r, km > V£- ?? = 0. 

Still with the same notation as above, we also define the normalized size of a section Sh(xo) as 

\L\ 2 

(2.15) a(<SUx )) = a(Q h (q )) :-- 



(detL) 2 /™' 

Notice that, even if L may not be unique, a is well defined up to universal constants. In case u 
is C 2 in a neighborhood of xq, by a simple Taylor expansion of u around qo it is easy to see that 
there exists h(xo) > small such that 

(2.16) a(S h (x )) = a(Q h (qo)) ~ \D 2 u(q )\ Vh < h{x ), 

where qo := q{xo). Since u and u are related by a diffeomorphism, the following lemma holds: 

Lemma 2.1. Let Q, 1 C and u £ C 2 (il / ) 6e a strictly c- convex function such that H-Dulli^^') is 
universally bounded. Then there exists a universal constants Mi such that the following holds: For 
every xq G $7 there exists a height h(xo) > such that if \D 2 u(xq)\ > Mi, then 

(2.17) \D 2 u(x )\ ^ a(S h (x )) Vh<h(x ). 
Proof. Differentiating twice the relation (|2.1|) we obtain 

D qq u = D q xD xx u + D qq xD x u + D q xD xx c + D qq xD x c, 

which implies that 

(2.18) u\D 2 u{q Q )\ - C(l + |Dn(xo)|) < \D 2 u{x )\ < ^-\D 2 u(q )\ + C(l + \Du(x )\) 

for some universal constants v, C > 0. Since by assumption D« is universally bounded inside f2', 
(12~T71) follows by flZHBD and (I27THD . provided Mi is sufficiently large. □ 
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We show now some geometric properties of sections and some estimates for solutions of (|1.4|) 
which will play a major role in the sequel. Here, the dilation of a section Sh(x) is intended with 
respect to x. 

Proposition 2.2 (Properties of section). Let u be a c-convex Aleksandrov solution of (ll.4p with 
< A < / < 1/A. Then, for any £1' CC Q" CC £1, there exists a positive constant p = p(Q',£l") 
such that the following properties hold: 

(i) S%(x) C n" for any x G SI', < h < 4p. 

(ii) There exist < ct\ < a.2 universal such that for all p G (0, 1) 

p^St{x) C Sl h (x) C p^Sl{x) 

for any x G Q' , < 2h < p. 
(hi) There exists a universal constant a < 1 such that, if S%(x) n S%(y) ^ 0, then S%(y) C 

S^i a {x) for any x, y G O', < h < up. 
(iv) n 0<h < p S%(x) = {x}. 

Proof. Points (i) and (iv) follow from the strict c-convexity of u shown in section 7] , and the fact 
that the modulus of strict c-convexity is universal (this last fact follows by a simple compactness 
argument in the spirit of [H Theorem 1']). 

Point (hi) corresponds the engulhng property of sections proved in [1 H Theorem 9.3]. 

The second inclusion in point (ii) follows from Lemma 9. 20. For the first one, it is enough 
to show that there exists a universal constant s£ (0, 1) such that 

(2.19) sQfM) C Qt /2 (q) 

and then iterate this estimate (here u is defined as in (|2.ip . and q := q(x)). To prove (|2.19p . let E 2 h 
be the John ellipsoid associated to Q^hi^^ an d assume without loss of generality that that E 2 h is 
centered at the origin. By convexity of the sections in this new variables, 

KQhti) QhiM) C QU<i) C nE 2h V s G (0, 1). 

Observe now that, for any q G Q^iq), we have (recall that u(q) = 0) 



(2.20) u(s(q - q) 



+ q)=s I Du{{\ -ts)q + tsq) ■ (q - q)dt. 
Jo 



Since q,q & nE 2 h we have q — q G 2nE 2 h, hence 

(2.21) {q-q)/2n€ E 2h cQ 2h {q). 

Moreover, by convexity of Q2h(<j), (1 — ts)q + tsq G Qh(<j) C ToQ 2 h(q) for some universal To < 1 
(see [HJ Lemma 9.2]). Defining the "dual norm" || • associated to a convex set K. as 

\\ a \\*K '■= supa • £, 

it follows from [1 H Lemma 6.3] that 

(2.22) \\Du(q)\\* Q2hm = || - D q c(q,T u (q))\\* Qk(q) < Ch \/ q G Qg(f). 



■*"To be precise, in [llj the dilation is done with respect to the center of the John ellipsoid, and not with respect 
to the "center" x of the section. However, it is easy to see that the same statement holds also in this case. 
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Thus, thanks to (gjZDJ) and $Z72T\i we get 

u(s(q -q)+q) = 2ns [ Du((l - ts)q + tsq) ■ — - - dt 
Jo 2n 

< 2ns J \\Du((l - ts)q + tsq))\\* Q2h ( q) dt < 2nsCh < h/2, 

provided s is small enough. This proves the desired inclusion. □ 

As shown for instance in [10] , an easy consequence of property (iii) is the following Vitali-type 
covering theorem. 

Proposition 2.3 (Vitali covering theorem). Let u, f, Q' , £1" , p, a be as in Proposition (12. 2p . let D 

be a compact subset offl', and let {Sh x (x)} x£ D be a family of sections with h x < p. Then we can 
find a finite number of these sections {Sh x .{ x i)}i=i,...,m such that 

m 

D C (J Sh x .(xi), with {Sah Xi ( x i)}i=i,...,m disjoint. 
i=i 

We now want to show that sections at the same height have a comparable shape. For this, we 
first recall the following estimate from 

Proposition 2.4 (Aleksandrov estimates). Let u, f,Cl' , O" , p be be as in Proposition ()2.2[) . and let 

Sh(x) be a section of u for some x E Cl' and h < p. Then 

(2.23) \S h {x )\^h n ' 2 . 

Remark 2.5. Estimates (|2.22p and (|2.23p have the following important consequence: consider 
the function u defined in (|2.ip . fix one of its sections Qh such that Qih C fi" with Jl" as above, 
normalize Qh using its corresponding John's transformation L, and define w as in (|2.10p . Since 
(detL)~ 2 / n w \Eh\ ~ oscQ h u w oscq 2Ii u (by (|2.23p ) and Eh C Q2h, we deduce the universal 
gradient bound 

sup \Dw\ = (det Lf' n supKL*)^ 1 Du\ 

L(Qh) Qh 

< C(detL) 2/n sup|| J Du||^ 

(2.24) Qh 

<C(detL) 2 / n S up\\Du\\* Q2h 

Qh 

< C(detL) 2 / n oscu < C. 



Lemma 2.6. Let u, f, £1' , O" , p be as in Proposition (|2.2p . Then for any < h < p there exist two 
radii r = r(h) and R = R(h) such that, for every xq E Q' , if E is the John ellipsoid associated to 
S^(xq), then, up to a translation, 

B r (0) Cfic5 fi (0). 

Proof Let r\ < . . . < r n be the axes of E. Since r n < diam(£') < C and by (I2.23P 

h n ' 2 «|£|»r r ...T„< diam( J E) n ^Vi, 
we obtain the desired lower bound on r\. □ 
Obviously analogous properties holds for the section Q^{qo)- 
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Remark 2.7. Notice that Proposition I2.2f ii) applied to the (convex) sections of u implies the 
following: given x G ft" and h < p, let r\ < . . . < r n denote the axes of the John ellipsoid 
associated to Qh(x). Then 

(2.25) r n < Crf 

for some universal exponent 03 < 1 and a constant C(£l',Q"). 

To see this just normalize Q p (x) using L and notice that, by [111 Theorem 6.11] , dist(x, d[L(Q p (x))) > 
1/C for some universal constant C. Thus, up to enlarge C, 

Bi(x) C L(Q h ) C f^M B 1 (x). 

Since, by Lemma 12.61 sections of height p have bounded eccentricity (i.e., \L\ ~ C(£l' this 
implies the claim with 03 := c^/ai • 

We now observe that cx{Qh) ~ ''n/ ( r i • • • r n) 2 ^ n , from which we deduce that 

r?<C- 5 



1_ a(Q h )' 
In particular, this and (|2.25p imply 

v' 

r n < Cr° 3 < C 



a(Q h )<*/*' 
that is 

r n <-^- R , with fi, ° 3 



<*W M " 2-2a 3 ' 

Hence, since is linked to Q/j by a diffeomorphism with universal C 1 norm, and diam(SVj) < 
diam(nEVi) = 2nr n , we get 

C 

(2.26) diamfS/,) < — 0, C > universal. 

3. H /2 ' 1+K ESTIMATES 

Applying first a large dilation to f2 we can assume that B(0, 1) C fi, and by a standard covering 
argument (see for instance [£l Section 3]) it suffices to prove the W 2,1+K regularity of u inside 
B(0, 1/2). Also, by an approximation argument^, it is enough to prove the result when u E C 2 . 
Hence Theorem II. II is a consequence of the following: 

Theorem 3.1. Let u G C 2 be a c-convex solution of (II. 4p with Q D 5(0,1). TTien t/iere e:ris£ 
universal constants k and C such that 

(3.1) I \D 2 u\ 1+K < C. 

JB(0,l/2) 

We start with the following lemma: 



9 

To approximate our solution with smooth ones, it suffices to regularize the data and then: 

- either apply [181 Remark 4.1] (notice that, by Proposition I2.2f iv) and Theorem 8.2], u is strictly c-convex and 
of class C 1 inside 

- or approximate our cost c with cost functions satisfying the strong version of (CI) and apply |18l Theorem 1.1]. 
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Lemma 3.2. Letu be as above, xq € 5(0,3/4), andh > such that S^/i^o) C 5(0,5/6). Consider 
the function u as in (12. ip . its section Qh = Qh(qo) with qo := q(xo), and (up to a rotation) 
let Eh = {Yl l x ll r f — l} ^ e ^ e </o/m ellipsoid associated to Qh(qo)- Denote by L be the affine 
transformation that normalizes Qh, and define w and Cij tk e as in (|2.10j) and ()2. 13j) respectively. 
Then 



(3.2) / dyit; - C ijtk £d k w d e w \< C 

JL{Q h ) 

for some universal constant C. 

Proof. Since, by the c-convexity of u (which is preserved under change of variables), the matrix 
(dijW — Cij kidkW dgw) . ._, is non-negative definite, it is enough to estimate 

„ n 

/ y2 i diiiB ~ Cii, s td k w diw) 

J L(Qh) ~1 

from above. 

Using the bounds V. 11 ' 1 (d[L(Qh))) < C(n) (since L(Qf l ) is a normalized convex set) and \Dw\ < 
C (see (|2.24p ). we see that first term is controlled from above by 

(3.3) I Avj= [ Dw-v<U n ~ x (d(L{Q h ))) sup \Dw\ < C. 

JL{Q h ) Jd(L(Q h )) L(Q h ) 

For the second term, we claim the following: there exists a universal constant C such that 



n 



(3.4) inf y^Cu.kedkVj dew > -C 
To see this we write 

n n n n 

C iiM d k w d e w = C iiM d k w d t w + 2 C iiy i k diW d k w + 2J Cu&diW diw. 

i=l i=l k,£y^i i=l k^=i i=l 

We first observe that, since for any i = 1, . . . , n the vector {d\w, . . . , di-\w, 0, c^+iif), . . . , <9 n w)) is 
orthogonal to coordinate vector ej, the first term in the right hand side is non- negative by condition 
(CI). 

Concerning the second and the third term, taking into account the definition of Cijki in (|2.12p 
we can rewrite them as 

n n 

(3.5) 2 ]T £( n . . . r n ) 2 ^B iitik diW 8 k w + . . . r n fl n B iiM diW 

i=l fc^i rfc i=l 

Observe that, by (f2~23l) . 

(3.6) (n...r n f/ fl »/i. 

In addition, by the Lipschitz regularity of u (which is simply a consequence of the fact that u is 
locally Lipschitz inside Jl), 

(3.7) h/r k <C Vfc = l,...,n. 

Since H-DtDH^ < C (see (|2.24|) ) and the size of £> is controlled by |||c||| ~ |||c|||, by (|3.6p and (|3.7p we 
see that the expression in (13. 5p is universally bounded. 

This proves (|3.4p . which combined with ()3.3p concludes the proof. □ 
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Lemma 3.3. With the same notation and hypotheses as in Lemma \3.''A let a be as in Proposition 
\2.3[ Then there exists a universal constant C such that 

(3.8) \{q G L{Q ah ) : U/C < d ijV d - C ijM d k wd e w < CId}\ > i. 

Proof. Since a is universal and L(Qh) is normalized, by Proposition I2.2f ii) we get 

\L{Qah)\ ~ \L(Q h )\ ~ 1. 

So, using Lemma 13.21 and Chebychev inequality, we deduce the existence of a universal constant C 
such that 

\{q G L(Q ah ) : dijW - C ijM d k w dgw < CId}| > — . 

Since by (|2.1ip the product of the eigenvalues of the matrix (pijW — Cij^dkw) i . =1 n is of order 
one, whenever the eigenvalues are universally bounded from above, they also have to be universally 
bounded also from below. Hence, up to enlarging the value of C, this proves (13. 8j) . □ 



Remark 3.4. Recalling the definition (12.150 of a(Qh) = cx.(Sh), we can rewrite both (|3.2|) and 
(13. 8p in terms of u and Qh = Qh(lo), obtaining that 



/ \diju - B ijM d k udiu\ < Ca(Q h ) {x G Q ah : a(Q h )/C < \diju - B ijM d k udiu\ < Ccx{Q h )} 
JQh 

(see for instance the proof of [101 Lemma 3.2]). In terms of u, this estimate becomes 

(3.9) / \D 2 u-A(x,Du)\ < C a(S h ) {x G S ah : a(S h )/C < \D 2 u - A(x, Du)\ < C a(S h )} , 

where 5^ = Sh(xo) with xq an arbitrary point inside 5(0,3/4), and Co is universal. 

Proof of Theorem \3.1[ Let Af ^ 1 to be fixed later, set Rq := 3/4, and for all m > 1 define 

(3.10) R m := i? m _i - CM-?. 

with C and /3 as in (|2.26p . Let us use denote B(0,R m ) by Bj^, set 

(3.11) F(x) := \D 2 u(x) - A{x,Du(x))\, 
and define 

(3.12) D m := {x G B Rm : F(x) > M m } . 

Thanks to Proposition 12.21 there exists p > universal such that Sh(x) C 5(0,5/6) for any 
x G 5(0,3/4) and h < 2p, and by Lemma 12.61 applied with h = p we get a(S p (x)) ~ 1. In 
addition, since (CO) implies that \A(x,Du)\ < C\ inside 5(0, 1) for some C\ universal, we see that 
\\D 2 u\ -F\ < d. Hence, using Lemma 12.11 we deduce that if M > M\ + C\ then there exists a 
small universal constant v > such that 

a(S h (x)) > vM m Vi G D rn , h < mm{h(x), p}, m > 1. 

So, by choosing M > max{l/i/ 4 , Mi} (so that vM m+1 > M m+ y 2 /v), by continuity we ob- 
tain that, for every point in D m+ i, there exists h x G (0, va.m.{h{x), p}) such that a(Sh x ) G 
(vM m+l / 2 ,M m+1 l 2 /u). In particular, by fl!T26l) we have diam^J < CM' 13 , which implies that 
(recall (131(1 ) 

(3.13) J S hx cB (o, 5 m+1 + CM-?) = B Rm . 

x£D m +l 
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According to Proposition 12.31 we can cover D m+ \ with finitely many sections {Sh x }xjeD m+1 such 
that S a h x . are disjoint. Then f|2. 1T|) and fj3.9j) imply (recall (|3.11j) ) 

/ F -Y,[ ^ — T> C ° a ^M x G S ° h * 3 : a ( S ^)/C <F< C oa (S hx .)}\ 

JD m+ i j JS hx , j 

< ^ — M m+l ' 2 \{x E S ahx . : uM m+1 ' 2 /C Q <F< C M m+1 / 2 /u}\. 

j 

Assuming now that VM > Cq/u and recalling ()3.13p . we obtain 
(3.14) 

/ F<Y [ F < V — M m+1 ' 2 \{x e S ah : vM m+1/2 /C <F< C M m+1/2 /u}\ 

< Y — M m+1 / 2 \{x eS ah :M m <F< M m+1 }\ 



9l M m+1 / 2 \{x E S ahx . :M m <F< M m+1 } D B Rv 



v 

j 



v JD m \D m+1 

Adding Q^M- J D ^ F to both sides of the previous inequality, we obtain 

(i+^H) ( F <9^Kl F . 

\ V ) JD m+1 v J Dm 

which implies 

F<(1-t) 



Dm+i J Dm 

for some small constant r = r(M) > 0. We finally fix M so that it also satisfies 

(3.15) CM~ ml3 < 

m>l 

In this way R m > 1/2 for all m > 1, so that the above inequalities and the definition of R m imply 



F< F<(l-r) m F<C(l-T) m 

{F>Af m }nB(0,l/2) J D m JDo 

(here we used that f B , Q 3 ^ F < C, which can be easily proved arguing as in the proof of Lemma 
3.2p . Thus, choosing k > such that 1 — r = M~ 2k , we deduce that 



\{F>t}nB(0,l/2)\<j [ 

1 h. 



F < Ct- l ~ 2K 

{F>t}nB(0,l/2) 



for some C > universal, which implies that F E L 1+K (B(0, 1/2)). Recalling the definition of F 
(see (15XT10 and that \A{x, Du)\ < C inside B(0, 1) (by (CO)), this concludes the proof. □ 



12 



G. DE PHILIPPIS AND A. FIGALLI 



References 

L. Ambrosio, M. Colombo, G. De Philippis, A. Figalli: Existence of Eulerian solutions to the semi- 
geostrophic equations in physical space: the 2- dimensional periodic case. Comm. Partial Differential Equations, 
to appear. 

L. Ambrosio, M. Colombo, G. De Philippis, A. Figalli: A global existence result for the semi-geostrophic 
equations in three dimensional convex domains. Preprint 2012. 

L.Caffarelli: A localization property of viscosity solutions to the Monge-Ampre equation and their strict con- 
vexity. Ann. of Math. (2), 131 (1990), no. 1, 129-134. 

L.Caffarelli: Some regularity properties of solutions to Monge-Ampere equations. Comm. Pure Appl. Math., 
44 (1991), 965-969. 

L. Caffarelli: The regularity of mappings with a convex potential. J. Amer. Math. Soc, 5 (1992), 99-104. 
M. J. P. Cullen, R. J. DOUGLAS, I. Roulstone, M. J. Sewell: Generalized semi-geostrophic theory on a 
sphere. J. Fluid Mech., 531 (2005), 123-157. 

P. Delanoe, Y. Ge: Regularity of optimal transport on compact, locally nearly spherical, manifolds. J. Reine 
Angew. Math., 646 (2010), 65-115. 

P. Delanoe, Y. Ge: Locally nearly spherical surfaces are almost-positively c-curved. Methods Appl. Anal., 18 
(2011), no. 3, 269-302. 

G. De Philippis, A. Figalli: W 2 ' 1 regularity for solutions of the Monge-Ampere equation. Invent. Math., to 
appear. 

G. De Philippis, A. Figalli, O. Savin: A note on interior W 2,1+E estimates for the Monge-Ampere equation. 
Math. Ann., to appear. 

A. Figalli, Y.H. Kim, R.J. McCann: Holder continuity and injectivity of optimal maps. Preprint, 2011. 
A. Figalli, G. Loeper: C 1 regularity of solutions of the Monge-Ampere equation for optimal transport in 
dimension two. Calc. Var. Partial Differential Equations, 35 (2009), no. 4, 537-550. 

A. Figalli, L. Rifford: Continuity of optimal transport maps and convexity of injectivity domains on small 
deformations o/§ 2 . Comm. Pure Appl. Math., 62 (2009), no. 12, 1670-1706. 

A. Figalli, L. Rifford, C. Villani: Nearly round spheres look convex. Amer. J. Math., 134 (2012), no. 1, 
109-139. 

F. John: Extremum problems with inequalities as subsidiary conditions. In Studies and Essays Presented to R. 
Courant on his 60th Birthday, January 8, 1948, pages 187-204. Interscience, New York, 1948. 
Y.-H. Kim, R. J. McCann: Continuity, curvature, and the general covariance of optimal transportation. J. 
Eur. Math. Soc. (JEMS), 12 (2010), no. 4, 1009-1040. 

J. Liu: Holder regularity of optimal mappings in optimal transportation. Calc. Var. Partial Differential Equations, 
34 (2009), no. 4, 435-451. 

J. Liu, N.S. Trudinger, X.-J. Wang: Interior C 2 ' a regularity for potential functions in optimal transportation. 
Comm. Partial Differential Equations, 35 (2010), 165-184. 

G. Loeper: On the regularity of solutions of optimal transportation problems. Acta Math., 202 (2009), no. 2, 
241-283. 

G. Loeper: Regularity of optimal maps on the sphere: The quadratic cost and the reflector antenna. Arch. 
Ration. Mech. Anal, 199 (2011), no. 1, 269-289. 

X. N. Ma, N. S. Trudinger, X. J. Wang: Regularity of potential functions of the optimal transportation 
problem. Arch. Ration. Mech. Anal., 177 (2005), no. 2, 151-183. 
T. Schmidt: W 2 ' 1+e estimates for the Monge-Ampere equation. Preprint 2012. 

N. S. Trudinger, X. J. Wang: On the second boundary value problem for Monge-Ampere type equations and 
optimal transportation. Ann. Sc. Norm. Super. Pisa CI. Sci. (5), 8 (2009) 143-174. 
[24] N. S. Trudinger, X. J. Wang: On strict convexity and continuous differentiability of potential functions in 
optimal transportation. Arch. Ration. Mech. Anal., 192 (2009), 403-418. 

SCUOLA NORMALE SUPERIORE, P.ZA DEI CAVALIERI 7, 1-56126 PlSA, ITALY 
E-mail address: guido.dephilippis@sns.it 

The University of Texas at Austin, Mathematics Dept. RLM 8.100, 2515 Speedway Stop C1200, 
Austin, Texas 78712-1202, USA 

E-mail address: figalli@math.utexas.edu 



